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I. INTRODUCTION 



Study of the strong Ec — > Acvr, S* ^ AcTt and S* — Scvr decays in a nonrelativistic 

quark model 

C. Albertus,'^ E. Hernandez,^ J. Nieves,^ and J. M. Verde- Velasco^ 

' Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, Spain. 
^Grupo de Fisica Nuclear, Facultad de Ciencias, E-37008 Salamanca, Spain. 

We present results for the strong widths corresponding to the Ec Ac n, E* — > Ac tt and H* — + 
Sc TT decays. The calculations have been done in a nonrelativistic constituent quark model with 
wave functions that take advantage of the constraints imposed by heavy quark symmetry. Partial 
conservation of axial current hypothesis allows us to determine the strong vertices from an analysis of 
the axial current matrix elements. Our results r(E++ A+ 7r+) = 2.41 ± 0.07 ± 0.02 MeV, r(E+ 
At 7r°) = 2.79 ± 0.08 ± 0.02 MeV, r{T,° A+tt") = 2.37 ± 0.07 ± 0.02 MeV, r(E*++ A+ 7r+) = 
. . . 17.52 ±0.74 ±0.12 MeV, r(E*+ ^A+7r°) = 17.31 ±0.73 ±0.12 MeV, r(E*° ^A+TT") = 16.90 ± 

in ' 0.71 ±0.12 MeV, r(H*+ ^ H°7r++H+7r") = 3.18 ±0.10± 0.01 MeV andr(H*" ^ H+tt^+H^tt") = 

' 3.03 ± 0.10 ± 0.01 MeV are in good agreement with experimental determinations. 

o ■ 

CN ■ PACS numbers: 11.40.Ha,12.39. Jh,13.30.Eg,14.20.Lq 

> 

o 

' The nonrelativistic constituent quark model (NRCQM), using QCD-inspired potentials, has proved to be an excel- 
, lent tool to predict properties of hadrons. In the case of baryons including one heavy quark c or 6 and two light ones 
J> ' M, d or s, we can take advantage of yet another property of QCD: Heavy quark symmetry (HQS) 0,0, 0,0- This 
^sD . symmetry arises when the heavy quark mass is much larger than the QCD scale (Aqcd)- In that limit the dynamics 
lO of the light quark degrees of freedom becomes independent of the heavy quark flavor and spin. The light degrees of 
freedom are thus well defined and the masses of the baryons depend only on the quark content and on the light-light 
■ quantum numbers. This simplification was used in Ref. p to solve the three-body problem for the ground state (L=0) 
jy-^ of baryons with a heavy quark using a simple variational ansazt. We obtained static properties (masses, mass and 
. charge radii. . . ) in good agreement with previous calculations that used more involved Faddeev equations Q. The 
' advantage of our approach is that we also obtain easy to handle wave functions. Those wave functions were already 
, used in Ref. to study, with good results, the semileptonic decays of Aft and Sf, baryons. 
II ' In this work we shall evaluate strong widths for the Sc Ac tt, E* ^ Ac tt and S* ^ Sc tt decays. Last decade has 
Qh! seen a great progress on charmed baryon physics and now the ground state baryons with a c quark, with the exception 
^ I of the JT! ) ^■re well established (8| , and we have experimental information on the strong one-pion decay widths for the 
; Sc Ullflllll, ESIIl and S* mill. With very little kinetic energy available in the final state these reactions 
^ 1 should be well described in a nonrelativistic approach. Although they have been analyzed before in the framework of 
the constituent quark model (CQM) [l5lll6|. no attempt was made there to evaluate the full matrix elements. While 
there have been dynamical calculations in other models (see references below), to our knowledge, ours is the first 
' dynamical calculation within a nonrelativistic approach. In our calculation we will use the HQS-constrained wave 
functions that we evaluated in Ref. j5] using different inter-quark interactions, and whose goodness have already been 
tested in the study of the semileptonic A;, — * Ac and — > Sc decays in Ref. ■ The use of different quark-quark 
potentials will allow us to obtain theoretical uncertainties on the widths due to the quark-quark interaction. The 
pion emission amplitude will be obtained in a spectator model (one-quark pion emission) with the use of partial 
conservation of axial current hypothesis (PCAC). 

These reactions, and similar ones, have also been addressed in QCD sum rules (QCDSR) p^.[T^. in heavy hadron 



chiral perturbation theory (HHCPT) 0, 0, |23, l2ll |2^ . and within relativistic quark models like the light-front 
quark model (LFQM) ^ and the relativistic three-quark model (RTQM) 



II. ONE-PION EMISSION AMPLITUDE ^^^^ (Ps , PsO 

To determine the pion emission amplitude A^^'^,1^{Pb, Pb') we shall use PCAC, as we have done in the meson 
sector, in a previous study of the strong B*Btt and D*Dit couplings [2^. PCAC allows us to relate that amplitude 
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to the matrix element of the divergence of the axial current. For the emission of a tt"*" we have-"^ 

{b\ s'pB'\q^jZiO)\B, sPb) ^ ^if.A^Bil+iPB.PB') (1) 

\ / non—pote 

where s, s' are the third component of the spin of the B, B' baryons in their respective center of mass systems, 
Pb — {.Eb{\Pb\), Pb), Pb' = {Eb'(\Pb'\), Pb') are their respective four-momenta and q ^ Pb - Pb'- "^1^(0) is the 
axial current for the u ~* d transition, and = 130.7 MeV Q is the pion decay constant. The baryon states are 

normalized as (^B, s'P' \B, sP^ = 5s. s' (27r)'^ 2Eb{\P |) S^{P — P'). Furthermore we shall use physical masses taken 
from Ref. in all calculations. 



III. DESCRIPTION OF BARYON STATES 



We use the following expression for the state of a baryon B with three-momentum P and spin projection s in the 
baryon center of mass 



1 



ai Pi 



Ql ,Q2 ,0^3 
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(27r)3 ^2Ef,m)2Ef,{\p2\)2Ef.Am) 
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M 
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"3 P3 



M 



P-Q1-Q2 



(2) 



aj represents the quantmn numbers of spin s, flavor / and color c (aj — {sjT fj,Cj)) of the j-th quark, while 
{Ef.(\pj\), pj) and rri/^. represent its four-momentum and mass. M stands for M = m/j + TO/2 + nif^. We choose 
the third quark to be the c quark while the first two will be the light ones. The normalization of the quark states 
IS { a' p ' \ a p) — 6a', a {2t^Y 2E{\p\) 5'^{p ' —p). Besides, 'i/'if 02,03 ( Qi, Q2 ) is the nonrelativistic momentum space 
wave function for the internal motion, being Qi and Q2 the momenta conjugate to the relative positions ri and f2 
of the two light quarks with respect to the heavy one. This wave function is antisymmetric under the simultaneous 
exchange ai < — > 02, Qi < — > Q2, being also antisymmetric under an overall exchange of the color degrees of freedom. 
It is normalized such that 



and, thus, the normalization of our nonrelativistic baryon states is 



NB. 



B,s' P'\B.sP 



NB. 



Ss>,s{27TfS^P' -P 



(3) 



(4) 



For the particular case of ground state Ac, Sc, S*, and S* we can assume the orbital angular momentum to be 
zero. We will also take advantage of HQS and assume the light-degrees of freedom quantum numbers are well defined 
(For quantum numbers see, for instance. Table 1 in Ref |^). In that case we have^ 



(Qi,Q2) 



(Qi,Q2) 



'(Ol,02) 



i^i^ (l/2,l/2,0;si,S2,0) 



V2 



{Sn,u €\tc^ Qi,Q2) -Sf,,d <5/2, n Ql, Q2 ) 



'C1C2C3 iSi = 



V3! 



ulciQu Q2 ) -5/,, „ 5/2^ „ <5/3,c 51(1/2, 1/2, 1; si, S2,m) (1, 1/2, 1/2; to, S3, s) 



-C1C2C3 iS,=l 



V3! 



'l^u:u!ciQi,Q2) Sf,,uSf,,uSf,,c ^(1/2, 1/2,1; Si, S2, to) (1,1/2, 3/2; to, S3, s) 



^ Note that we give the expression corresponding to the non— pole part of the matrix element. If the pion pole contribution is included 

then the relation is given by (S', s' Pg' \ qt" Jl^(O) | B, sPg ) = -if-^-^^ A^^'^^l^iPg , Pb')- 
^ We only give the wave function for the baryons involved in 7r+ emission. Wave functions for other isospin states of the same baryons 

are easily constructed. 
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U:hto.s{Qi,Q2) = (1/2, 1/2,0; si,S2,0) 

X ^-^^^^ [Sf,^, 6f,,s 4>l7'AQi,Q2 )~Sf,,,Sf,^, i%:°iQi,Q2 )) 

4'&!UQi,Q2) - ^(1/2, 1/2,1; SI, S2,m) (1,1/2, 3/2; m, S3, s) 

X ^ fe^(Qi,Q2) + <5/„.%,„ fc'(Qi,Q2)) (5) 

Here £ciC2C3 is the fully antisymmetric tensor on color indices being 6010203/ Vsl the antisymmetric color wave function, 
the (j'l, j2, j; TOi, '712, ma) are Clebsch-Gordan coefficients and the (f)j^^ f^QiiQi), with Si the total spin of the 
light degrees of freedom, are the Fourier transform of the corresponding normalized coordinate space wave functions 
obtained in Ref. j^. Their dependence on momenta is through |Qi|, IQ2I and Qi ■ Q2 alone, and they are symmetric 
under the simultaneous exchange /i < — > 12, Qi < — ^ Q2- 



IV. RESULTS 

Ec — * Ac TT decay 



Let us do the A+ 7r+ case. The matrix element of the divergence q'^J'A^iO) of the axial current determines 

the 7r+ emission amplitude as 



A'^^:l{^^^iP,p') ^ ^ (a+, s'p' |9^4;i(o)| sp) ^^g^t^At.+ ^^t Ap')i^ ^e^+ .(^) (6) 

where the coupling constant 5s++A+7r+' analogy to the pion coupling to nucleons and nucleon resonances, has been 

chosen to be dimensionless, and g{P), Up^+ gt{P') are Dirac spinors normalized to twice the energy. The width 

is given by 

r(S++ ^ A+.+ ) = al,,^,^, ( (M^.+ - M^,f - ml) (7) 

with \q \ the modulus of the final baryon or pion three-momentum. From Eq. lO, taking P ~ 0, P' ~ —\q \ k in the z 
direction, s = s' = 1/2, and taking into account the different normalization of our nonrelativistic states, we have 



5s++A+7r+ = -r 



with 



((Af,.. - E,^m) 4/^4/^^ ^ + 1,-1 A^/^^/J^ 3) (8) 



4t+i+,;. = A'«(A- 1/2 -\q]k\ JiliO) I 1/2 0)..^ , (9) 



NR, non-pole 



The Aj^^ are easily evaluated using one-body current operators and their expressions can be found in the 



,1/2,1/2 

S++A+, p 

appendix. 

In Table n we present the results for g^++j^+^+ and the widths r(S++ A+7r+), r(S+ A+7r°) and r(I]^ -> 
A+TT^). To get the values for r(S+ A+tt") and r(E" — > A+vr^) we use 35]++A+7r+ ^^'^ make the appropriate mass 
changes in the rest of factors in Eq. |7J. Our results show two types of errors. The second one results from the Monte 
Carlo evaluation of the integrals needed to obtain the .95]++A+7r+ coupling constant. The first one, can be considered as 
a theoretical uncertainty and comes from the use of different quark-quark interactions ^. The results are in very good 



^ We use five different int er-quark interactions, one suggested by Bhaduri and collaborators l2fill . and four otfiers suggested by Silvestre- 
Brac and Semay I^ l27l . All of them contain a confinement term, plus Coulomb and hyperfine terms coming from one-gluon exchange, 
and differ from one another in the form factors used for the hyperfine terms, the power of the confining term or the use of a form factor 
in the one gluon exchange Coulomb potential. 
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This work 



Experiment 



Theory 
CQM 

HHCPT 



LFQM 
RTQM 



" [MeV] [MeV] [MeV] 

21.73 ±0.32 ±0.08 2.41 ± 0.07 ± 0.02 2.79 ± 0.08 ± 0.02 2.37 ± 0.07 ± 0.02 



2.3 ± 0.2 ± 0.3 [9] < 4.6 (CL=90%) [lOj 2.5 ± 0.2 ± 0.3 [9J 
2.05lg;38 ± 0-38 Hi] 1.55lg;37 ± 0.38 [llj 



1.31 ± 0.04 [15] 1.31 ± 0.04 [15] 1.31 ± 0.04 [15j 

2.025+Jir, [16] 1.939^JiJl [16] 

22, 29.3 [19] 2.47, 4.38 [19] 2.85, 5.06 [19] 2.45, 4.35 [19] 

2.5 [22] 3.2 [22] 2.4 [2^] 

1.94 ±0.57 [21] 
1.64 [2^ 1.70 [2^] 1.57 [2^ 
2.85 ± 0.19 [24] 3.63 ± 0.27 [24] 2.65 ± 0.19 [24] 



TABLE I: Couphng constant g^++j^+^_^_ and total widths r(E++ A+7r+), r(E+ A+vr") and r{T^° A+vr") (See text for 
details). Experimental data and different theoretical calculations are also shown. 

agreement with the experimental data by the CLEO Collaboration in Refs. [ ft The value for r(I]++ A+7r+) 
also agrees with the experimental data by the FOCUS Collaboration in Ref. [ij- The agreement with FOCUS data 
is not good for the r(5]° A+7r~) case, although our result is still within experimental errors. Our results show 
variations as large as «17% between different charge configurations. This is due to the little kinetic energy available 
in the final state that makes the widths very sensitive to the precise masses of the hadrons involved. In this respect 
there is a new precise determination of the A+ mass by the BABAR collaboration Ma^ = 2286.46±0.14 MeV/c^ [28| . 
which is roughly 1.5 MeV/c^ above the value quoted by the Particle Data Group in Ref. With this new value 
our calculated widths would get reduced by 9%. This reduction comes from phase space factors while the coupling 
g^++^+^^ changes only at the level of 0.1%. 

As for the other theoretical determinations, the CQM calculation in Ref. uses exact s < — > c symmetry to relate 
the Sc — > AcTT decay to the non-charmed S* — > Att analogue decay. Their results are smaller than the experimental 
data obtained by the CLEO Collaboration. In Ref. 16] a unique coupling constant is fixed in order to reproduce all 
experimental information on E* AcTt widths. That coupling is latter used to predict the Ec ^ Acvr widths. This 
coupling suffers from large uncertainties and thus the theoretical errors on the predicted widths are also very large. In 
the HHCPT calculation of Ref. 19] a simple CQM argument is used in order to obtain the unknown coupling constant 
in the HHCPT Lagrangian. Furthermore the authors allow for a renormalization of the axial coupling g'^ for light 
quarks. The largest of the two values quoted corresponds to the case where that coupling is unrenormalized and then 
is given by = 1. The smaller number quoted corresponds to the use of a renormalized value of g^ — 0.75. The 
case gj^ = 1 is the one that compares with our calculation. Their results for the widths almost double ours and are not 
in agreement with experiment. Their simple determination of the coupling constant can not be correct. The values 
obtained in the HHCPT calculation of Ref. [23] are closer to our results and experimental data. There the authors 
determine the needed coupling constants from the analysis of analogue decays involving non-charmed baryons. In 
Ref. |21|, also within the HHCPT approach, and similarly to Ref. jl^, the authors fix the unknown coupling in the 
Lagrangian using the experimental data for the S* AcTt decays. From there they predict the r(I]" A+tt^) 
obtaining a value very close to the one in Ref. 16], and that suffers also from large uncertainties. The two relativistic 
quark model calculations of Refs. 23, 24] give results that differ by almost a factor of two. Our results are closer to 
the ones obtained within the RTQM of Ref. 

B. E* — > AcTT decay 

Let us analyze the case with a 7r+ in the final state, A+7r+. Similarly to the Sc decay before we now have 

4r;i^,^^(P,P') = ^(a+ s'P' |g^4^(0)| sP) 
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r(Sc AJ TT ) 








[MeV] 


[MeV] 


[MeV] 


This work 


36.20 ±0.75 ±0.13 


17.52 ±0.74 ±0.12 


17.31 ±0.73 ±0.12 


16.90 ± 0.71 ± 0.12 


Experiment 




14.1+i;| ± 1.4 [12] 


< 17 (CL=90%) [10] 


IG.&ti'j ± 1-4 [12] 


Theory 










QCDSR 


13.8-^24.2 [17] 










32.5 ±2.1 ±6.9 [18] 








CQM 


20 [15] 


20 [15] 


20 [15] 


HHCPT 




25 [20] 


25 [20] 


25 [20] 


LFQM 




12.84 [2^ 




12.40 [2^ 


RTQM 




21.99 ± 0.87 [24] 




21.21 ±0.81 [24] 



TABLE II: Coupling constant g^,+^^+^^ and total widths r(E*++ A+7r+), r(E*+ A+7r°) and r(E*° ^ A+vr"). 
Experimental data and different theoretical calculations are also shown. Note that in order to compare with our definition of 
g^, we have multiplied the coupling constants evaluated in Refs. fr^ ll8|l by 2M^+ //^. 

where we have introduced the dimensionless coupling constant 5s*++a+^+ and ^{P) is a Rarita-Schwinger 

spinor normalized to twice the energy. The width is given by 

2 

\q? gs;++A+7 r+ 

Taking again P = Q, P' ~ —\q\k in the z direction, and s = s' = 1/2 we obtain from Ea. (|10|l 



5s*+ + A+7r+ 



with 



^3 2M^^J2M^,^+2E^^{\q\) 



UV2 



q\j2M^.^+{E^U\q\) + M^ 



((Afj,.+. ~ E^,{\q\)) 4(';V1+, + k1 4^';VI+, 3) (12) 



41/2,1/2 

e;++A+, n NR 



The expressions for A^^^lY^^. (/i = 0, 3) can be found in the appendix. 
Results for 5^^. ++^+^+ and the total widths r(S*++ A+7r+), r(E*+ A+tt") and T{J:*- A+tt") appear in 
Table HTl Our value for the latter two are obtained with the use of 5x;*++A+Tr+ ^n*^ with the appropriate mass changes 
in the rest of factors in Eq. (|ll|l . Our central value for r(E*"'"+ A+tt"*") is above the central value of the latest 
experimental determination by the CLEO Collaboration in Ref. For some of the potentials used, API and AP2 
of Ref. 0, the results obtained are within experimental errors. The central value for r(E*+ A+7r°) is slightly 
above the upper experimental bound determined also by the CLEO Collaboration in Ref. 10], but again, we obtain 
results which are below the experimental bound using the API and AP2 potentials. As for r(E* ^ A+tt^) wc get 
a nice agreement with experiment. Our results for the different charge configurations differ by 4% at most. With the 
new value for given by the BABAR Collaboration in Ref. l23| they would get reduced by 3%. Our results are 

globally in better agreement with experiment than the ones obtained by other theoretical calculations with perhaps 
the exception of the QCDSR calculation of Ref. 0. 



(a+, 1/2 I 4;(0) I 1/20) (13) 

\ / N non—pote 



C. He ScTT decay 

Once more we analyze the case with a 7r+ in the final state, S*+ ^ ^c''^^ ■ What we obtain is 



A'::iL^4P,p') = -(so, s'p'\q'^Jii{o)\E:+, sp 



^ We did not show those cases were data on S* AcTv widths were used to fit parameters of the models. 





[MeV] [MeV] 


[MeV] [MeV] 


This work 


-28.83 ±0.50 ±0.10 1.84 ± 0.06 ± 0.01 1.34 ± 0.04 ± 0.01 


2.07 ±0.07 ± 0.01 0.956 ± 0.030 ± 0.007 


Theory 
LFQM 
RTQM 


1.12 [23| 0.69 [23| 
1.78 ± 0.33 [Ml 1.26 ± 0.17 [24] 


1.16 [23] 0.72 [23] 
2.11 ± 0.29 [24] 1.01 ± 0.15 [24] 




r(H:+^H«^+±H^7r°) 
[McV] 


r(H*+-H^+±H,+^°) 
[McV] 


This work 


3.18 ±0.10 ±0.01 


3.03 ±0.10 ±0.01 


Experiment 


< 3.1 (CL=90%) fl3] 


< 5.5 (CL=90%) [14] 


Theory 
CQM 

HHCPT 

LFQM 

RTQM 


< 2.3 ±0.1 [15] 
1.191- 3.971 [1£| 
2.44 ± 0.85 [21] 

1.81 [23] 
3.04 ± 0.50 [24] 


< 2.3 ±0.1 |15] 
1.230 - 4.074 [1^ 
2.51 ±0.88 [21] 

1.88 [2^ 
3.12 ± 0.33 [24] 



TABLE III: Values for the coupling 5=.+20^+ and decay widths r(H*+ H^7r+), r(H*+ -* E+tv°), r(H*° H+tt") and 

r(H*'' — > HJ^tt''). Experimental upper bounds for the total H*+ and H*° widths, and different theoretical calculations are also 
shown. 



s^^-^^; 9. %o AP') (14) 



+ a/ho 

where again we have introduced a dimensionless coupling g^t+^g^^. The width is given as 



Taking now P = 0, P' — —\q \ k in the z direction, and s = s' = 1/2, (7„.+„o^+ is evaluated from Ea. H14|l to be 
^ (A%+MHo)y2Af...2i?Ho(|g|) 



UV2 \q\ /2M^,+ (E^o{\q\) + M^o] 



. = nr(^c, 1/2 - i9ifci 4';i(o) 1/20),,^ , a?) 



with 

41/2,1/2 

which expressions can be found in the appendix. 

Results for the coupling gH-+HO^+, the widths r(S*+ ^ S07r+), r(S*+ ^ S+tt"), r(S*o -> S+tt") and r(S*o ^ 

S^tt"), and the total widths T{k*+ E°tt+ + S+7r°) and r(S*o ^ E^^tt^ + S+tt") appear in Table UTTI Our values 
for r(S*+ ^ S+ttO), r(S*o S+tt") andr(S*° ^ SOtt") are obtained with the use of g^.+^o^+ /{M^+ + M^o), and 

with the appropriate mass changes in the rest of factors in Eq. lfT5|l . For r(S* + 2+7r°) and r(S* ° S°7r°) an extra 
1/2 isospin factor should be included. Our results for r(S*+ S"7r+ + S+7r°) are slightly above the experimental 
bound obtained by the CLEO Collaboration As for the E* decay case above, the API and AP2 potentials 

gives results closer to experiment. For r(S*° — > S+vr" ± S"7r°) our result is well below the CLEO Collaboration 
experimental bound in Ref. (14| . Isospin breaking due to mass effects is clearly seen when comparing the predictions 
for r(S*+ S^vr^) and r(S*+ E!+7r°). One finds a factor 1.4 difference when a factor of two would be expected 
from isospin symmetry. Again, the fact that there is little phase space available makes the results very sensitive to 
the actual mass values. Compared to other theoretical calculations our results agree nicely with the ones obtained 
within the RTQM of Ref. while they are larger than most other determinations. 
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V. CONCLUDING REMARKS 



We have evaluated the widths for the charmed-baryon decays S* — > AcTt, Sc — > AcTt and S* — > ScTt within the 
framework of a nonrclativistic quark model. While there have been dynamical calculations of these reactions in 
other models, to our knowledge this is the first time that such a calculation has been done within a nonrclativistic 
approach. We have used wave functions constrained by HQS and that were obtained solving the nonrelativistic three- 
body problem with the help of a simple variational ansazt. For that purpose we took five different nonrclativistic 
quark-quark interactions that included a confining term plus Coulomb and hyperfine terms coming from onc-gluon 
exchange. To evaluate the pion emission amplitude we have used a spectator or one-quark pion emission model. The 
amplitude has been obtained with the use of PCAC from the analysis of weak current matrix elements. Our results 
are rather stable against the quark-quark interaction, with variations in the decay widths at the level of 6-8%. We 
find an overall good agreement with experiment for the three reactions. This agreement is, in most cases, better than 
the one obtained by other models. 
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APPENDIX A: EXPRESSIONS FOR THE J^l'i:^''^ MATRIX ELEMENTS 



The values for the ^^b^^^ evaluated using one-body current operators and their expressions are given by 

X 5:(l/2,l/2,l;.„-s„0)(l/2,l/2,0;.„-.„0) ^AQr ' W\k) 1,1. u^sAQ.) ^^1) 

^j2Ed{\Qi-\q\k\)2Eu{\QA) 

where the quark Dirac spinors are normalized to twice the energy. For = 0, 3 we get the final expressions 



(Ed{\Qi-\q\k\)+m^ + m„) 

2Ei{\Qt-\q\k\)2Eu{\Qi\) 



Ql 



Qi - m 



Eu{\Qi\) + mu Ea{\Qi-\q\k\) + md 



q\k, Q2 + = — \q\k) 



Ed{\Qi-\q\k\)+md^ (^Eu{\Qi\)+mu) / ^ 2{Qlf - \Qi\^ - Ql\q 



2Ed{\Qi-\q\k\)2E^{\Qi\) 



(^Ed{\Qi - \q\ k\) + rridj |) + m„) 



(A3) 



For A^{^4.y^+ ^ we just have 



,1/2,1/2 



-V2A 



1/2,1/2 



(A4) 



Similarly we get for AI/^iI(^ , = 0, 3 



41/2,1/2 

^H!+HO 



^ I d^gi d^Q^ ^f;=^(Qi, Q2) [0tT°(^3i - ^^^=;^ k1 fc, Q2 + =^ |g I fc) 



(^Ea{\Qi~\q\k\)+md) (s„(|Qi|) 
2i;d(|Oi-|g|fc|)2£;„(|Qi|) 



Ql 



+ 



Ql - 1? I 



,K(|Qi|) + m„ £;d(|gi - |g|A;|) + md^ 



91^, Q2 + =^\q\k) 



(Ea{\Qi-\q\k\)+md) (Eu{\Qi\)+mu) 
2Ei(\Q^-\q\k\)2E^{\Qr\) 



2m)-\Qi\-Qim 

{Ed{\Qi-\q\k\)+m^ (K(|Qi|)+to«) 
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